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Heteroduplexes formed between 
different members of a gene family 
migrate more slowly in non-denaturing 
polyacrylamide gels (PAGE) than 
bomoduplexes. Heteroduplexes with 
deletions/insertions as small as a 
single nucleotide are identified by 
a mobility retardation in PAGE. 
Heteroduplexes containing only 
mismatches also display a mobility shift . 
whose magnitude* depends on the % 
sequence divergence. Mobility shifts 
of deletion/insertion heteroduplexes 
are affected by the size, location and 
number of ddetkms/insertions and 
by the sequence of the unpaired and 
neighboring mismatched nucleotides. 
The heteroduplex mobility assay (HMA) 
and heteroduplex tracking assay (HTA) 
of the present invention can be used 
for tracking viral quasispecies evolution 
within and between individuals, for 
tracking species of infectious agents, and 
for monitoring DNA sequence changes 
including identification and tracking of 
cellular DNA polymorphisms. 
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A HETERODUPLEX MOBILITY ASSAY FOR THE ANALYSIS 
OP NUCLEIC ACID SEQUENCE DIVERSITY 

Field of the Invention 

5 frhe present invention describes a method for 

genetic analysis of gene families: the 
Heteroduplex Mobility Assay (HMA) . The method has 
utility in the analysis of gene pools, in 
particular, gene pools of disease causing 

10 microorganisms . Further, the present invention 

describes a method for evaluating the effects of a 
disease treatment protocol on DNA sequence 
variation of a nucleic acid target sequence 
associated with the disease, using tKe 

15 Heteroduplex Mobility Assay (HMA) or Heteroduplex 
Tracking Assay (HTA) . The method has utility in 
evaluating sequence variation among ^erie pools of 
disease causing microorganisms and variant 
disease-related genes, such as oncogenes • 
20 : . ! " 
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Background of the invention 

20 The ability to detect small changes in DNA 

sequences, such as, base substitutions, deletions, 
and insertions, has become important in a number 
of applications including diagnosis of human 
genetic diseases. Several methods for the 

25 detection of small sequence changes have been 
proposed. 

One method involves the use of restriction 
fragment polymorphism (RFLP) analysis. 
Restriction fragment polymorphisms arise when a 

30 base change results in the loss or acquisition of 
a restriction endonuclease cleavage site in a 
defined DNA sample, for example, a selected region 
of a viral genome. However, since the chances are 
low that a particular base change will eliminate 

35 or create a restriction site, RFLP analysis is 
labor intensive and tends not to provide a high 
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level of sensitivity since most changes in a DNA 
sequence go undetected. 

Two methods for the detection and 
localization of base changes have been described 
by .Myers, etal. (1987). In both methods, a 
radioactively labeled single-strand wild-type 
sequence probe is annealed to a test sample of 
cloned or genomic DNA. if the sample DNA contains 
an base change relative to the wild-type sequence, 
then a mismatch is formed at that site, in the 
first method, a DNA or RNA probe is used and 
mismatched duplexes are separated from perfectly 
paired duplexes by denaturing gradient gel 
electrophoresis (Myers, et al., 1985a), In the 
second method, the probe is a single-strand RNA 
molecule having a defined sequence. In this 
method ribonuclease is used to cleave any 
mismatched duplexes. Any cleaved products are 
then identified by polyacrylamide gel 
electrophoresis and autoradiography (Myers, et 
al., 1985b). 

Polymerase chain reaction has provided the 
means to obtain large scale amplification of 
target DNA (Mullis; Mullis, et al.) from, for 
example, defined genomic regions of DNA. A number 
of approaches have been utilized, for the 
detection of DNA mutations, in combination with 
polymerase chain reaction amplified products. 
These approaches include denaturing polyacrylamide 
gel electrophoresis (see above and Sheffield, et 
al., 1989), single-chain conformation polymorphism 
analysis (Orita, et al.), chemical cleavage of 
mismatches (Cotton, etal.). Constant denaturant 
gel electrophoresis has also been employed as a 
screening method for mutations in a defined DNA 
sequence (Borresen, et al., 1991) 



summary of the Invention 

The present invention provides a method of 
evaluating sequence diversity in a mixture of 
nucleic acids containing a target sequence. In 
the method of the present invention, amplification 
primers are selected which are complementary to 
nucleic acid sequences flanking the target region. 
The nucleic acids and the primers are combined 
under conditions that promote the hybridization of 
the primers to the nucleic -acids, thus generating 
primer/ nucleic acid complexes. These complexes 
are converted to double-strand fragments in the 
presence of a suitable polymerase and all four 
deoxyribonucleotides. The primer-containing 
fragments are amplified by repeated rounds of 
primer extension until a desired degree of 
amplification has been achieved. The resulting 
double-stranded amplified fragments are denatured 
and renatured to form a population of amplified 
fragment DNA duplexes. These duplexes are 
separated on polyacrylamide gels. The relative 
migration of the duplexes is analyzed to establish 
the relative degree of sequence relatedness in the 
population of amplified fragments. 

In the method, the starting nucleic acids can 
be RNA molecules which are converted to DNA 
templates using reverse transcriptase. 

Typically, the denaturing step of the method 
is thermal denaturing, and double-stranded 
fragments are generated using a thermostable DNA 
polymerase. 

In one embodiment, the amplified duplexes are 
separated by polyacrylamide gel electrophoresis 
(PAGE), and analyzed by visualization of the 
amplification products with ethidium bromide 
staining or autoradiography. 



In another embodiment, the amplified duplexes 
are separated by PAGE. The locations of the 
duplexes are then analyzed by transferring the 
nucleic acids from the gel to a support membrane, 
and hybridizing the nucleic acid transferred to 
the ^membrane with a labelled probe specific for 
the desired amplification products. 

In the present invention, the primers may 
contain at least one detection moiety* The 
detection moiety can be, but is not limited to, 
one or more of the following: a radioactive 
moiety, biotin, digoxigenin, or a chemiluminescent 
moiety. 

7 The method of the present invention can be 
applied to the analysis of nucleic acid from any 
infectious agent, including viruses> bacteria, 
mycoplasma, parasites, and fungi. One exemplary 
application of the present method is in evaluating 
intra- and inter-patient sequence diversity of 
Human Immunodeficiency Virus (HIV) . Further, the 
present method can be applied to the analysis of 
sequence variation between, def ined genetic loci, 
such as, oncogenes , protooncogenes , and disease- 
associated loci (e.g., Duchenne's muscular 
dystrophy) . This method^ is also useful to 
evaluate sequence diversity over time in mixtures 
of nucleic acids containing a target sequence, 
where the mixtures are serially obtained from a 
single source, for example, the single source may 
be a patient infected with HIV-1* 

The present invention also includes a method 
of evaluating sequence diversity between two 
different sample mixtures of nucleic acids, where 
the nucleic acids contain a target sequence. In 
this method, the two samples are treated as 
described above. After amplification, the two 
amplified samples are mixed, denatured, and 
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renatured. The resulting duplex molecules are 
separated on pblyacrylamide gels. The relative 
migration of the duplexes is then analyzed to 
establish the relative degree of sequence 
5 relatedness among the amplified fragments of the 
population. This method can be applied to the 
analysis of samples from different geographic 
locations, different patients, or for the same 
patient with different samples collected over 
10 time. 

In one embodiment of this method, the 
amplified fragments of one sample nucleic acid are 
labelled with a detection moiety and the labeled 
fragments are mixed with a molar excess of the 

15 amplified fragments of the other sample. 

Yet another embodiment of the present 
invention is a method for detecting the presence 
of a selected nucleic acid target region in a 
nucleic acid sample. In this method a duplex DNA 

20 probe having two complementary strands is 

selected, where the duplex is homologous to the 
target region arid each strand contains a detection 
moiety. Amplification primers complementary to 
nucleic acid sequences flanking the target region 

25 of the nucleic acid are also selected. The target 
region is amplified. The amplified products are 
denatured and mixed, in molar excess, with the 
duplex probe. The mixture is denatured, 
renatured, and analyzed as described above to 

30 establish the relative degree of sequence 

relatedness between the probe and sample target 
regions. 

Typically, the mixing ratio for amplified 
fragments to probe is 100:1, amplified fragments 
35 to probe. 

The method of the present invention has been 
described with reference to using amplification to 
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obtain samples suitable for heteroduplex mobility 
•analysis. However, other methods of obtaining 
nucleic acid samples can be used as well, 
including, but not limited to, cloning sequences 
of interest or isolating particular restriction 
endonuclease digestion products from genomic DNA 
preparations. 

Further, the present invention prbvides a 
method of evaluating, in mixtures of nucleic 
acids, the effect over time of a disease > j 
treatment, on DNA sequence variation of a nucleic 
acid target sequence associated with the disease. 
In the method of the present invention, 
amplification primers are chosen that are c 
complementary to nucleic acid sequences flanking 
the target sequence. Mixtures of nucleic acids 
are serially obtained from a single source. For 
example, sera is collected from an HIV-infected 
individual before treatment and at selected time 
points over the course of treatment. 

The nucleic acids from each sample and the 
primers are combined under conditions that promote 
the hybridization of the primers to the nucleic 
acids, thus generating primer/nucleic acid 
25 complexes. For each sample, these complexes are 
converted to double-strand fragments in the 
presence of a suitable polymerase and all four 
deoxyribonucleotides. The primer-containing 
fragments are amplified by repeated rounds of 
30 primer extension until a desired degree of 
amplification has been achieved. 

The resulting double-stranded amplified 
fragments are denatured and renatured to form a 
population of amplified fragment DNA duplexes for 
35 each sample. These migration of the duplexes on 
the gel is analyzed for each sample to determine 
the relative migration of the duplexes and 
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establish the relative degree of sequence 
relatedness in the population of amplified 
fragments. The effect of the treatment is 
evaluated by comparing the relative degree of 
5 sequence relatedness of amplified fragments in 
each serial sample between the serial samples. 
Particularly, comparing to the pre-treatment, or 
zero time point. 

In the method, the starting nucleic acids can 

10 be RNA molecules which are converted to DNA 
templates using reverse transcriptase. 

Typically, the denaturing step of the method 
is thermal denaturing, and double-stranded 
fragments are generated using a thermostable DNA 

15 polymerase. 

In one embodiment, the amplified duplexes are 
separated by polyacrylamide gel electrophoresis 
(PAGE) , and analyzed by visualization of the 
amplification products with ethidium bromide 

20 staining or autoradiography. 

In another embodiment, the amplified duplexes 
are separated by PAGE. The locations of the 
duplexes are then analyzed by transferring the 
nucleic acids from the gel to a support membrane, 

25 and hybridizing the nucleic acid transferred to 
the membrane with a labelled probe specific for 
the desired amplification products. 

In the present invention, the primers may 
contain at least one detection moiety. The 

30 detection moiety can be, but is not limited to, 
one or more of the following: a radioactive 
moiety, biotin, digoxigenin, or a chemi luminescent 
moiety. 

The method of the present invention can be 
35 applied to the analysis of nucleic acid from any 
microorganism/ infectious agent, including viruses, 
bacteria, mycoplasma, mycobacteria, parasites, and 
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fungi, one exemplary application of the present 
method is in evaluating intra- patient sequence 
diversity of Human Immunodeficiency virus (HIV) . 

Further, the present method can be applied to 
the analysis of sequence variation between defined 
genetic loci, such as, oncogenes, protooncogenes , 
and disease-associated loci {e.g., Duchenne's 
muscular dystrophy) . 

Yet another embodiment of the present 
invention is a method for evaluating the effect of 
a disease treatment procedure on the presence of a 
selected nucleic acid target region in a nucleic 
acid sample, in this method a duplex DNA probe 
having. two complementary strands. is selected, 
where the duplex is homologous to the target 
region and each strand contains a detection 
moiety. 

Serial samples are obtained from a single 
source, e.g., a first sample before treatment and 

20 a second sample after treatment. Amplification 
primers complementary to nucleic acid sequences 
flanking the target region of the nucleic acid are 
also selected. The target region is amplified in 
each sample. The amplified products from each 

25 sample are denatured and mixed, in molar excess, 
with the duplex probe. Each mixture is denatured, 
renatured, and analyzed as described above to 
establish the relative degree of sequence 
relatedness between the probe and sample target 

3 0 regions . 

Typically, the mixing ratio for amplified 
fragments to probe is 100:1, amplified fragments 
to probe. 

The method of the present invention has been 
described with reference to using amplification to 
obtain samples suitable for heteroduplex mobility 
analysis. However, other methods of obtaining 
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nucleic acid samples can be used as well, 
including, but not limited to, cloning sequences 
of interest or isolating particular restriction 
endonuclease digestion products from genomic DNA 
5 preparations. 

These and other objects and features of the 
invention will be more fully appreciated when the 
following detailed description of the invention is 
read in conjunction with the accompanying examples 
10 and drawings. 

Brief Description of the Figures 

Figure IB shows the result of heteroduplex 
mobility (HIV-l envelope gene DNA) analyzed on a 

15 2.5% agarose gel. Figure 1A shows the same , 
heteroduplex reactions separated on a 5% 
polyacrylamide gel. The heteroduplex samples were 
obtained from nested PCR amplifications performed 
using uncultured HIV-l seropositive PBMC DNA. 

20 Figures 2A to 2E show the results of 

heteroduplex analysis carried out with gap- and 
base-pair-mismatch-containing heteroduplex 
molecules. The heteroduplexes were resolved on 
agarose (2A) and polyacrylamide (2B-2E) gels. 

25 Figure 3A shows a plot of the mobility shifts 

of heteroduplexes against the % divergence , based 
on nucleotide sequence analysis, between the two 
strands of the heteroduplexes. Figure 3 B shows 
some of the heteroduplex mobility assay data used 

30 to generate the plot presented in Figure 3A. 

Figure 4 shows the results of a heteroduplex 
mobility assay analysis of variant sequences from 
different HIV-l quasispecies. In the figure, MA 
is HIV-positive asymptomatic, PE and BU were HIV- 

35 positive and symptomatic. 

Figure 5 shows polyacrylamide heteroduplex 
mobility of inter-individual DNA heteroduplexes 
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using complex DNA mixtures : (quasispecies) as probe 
and target, in the figure, "*" is unreannealed 
single stranded DNA. US quasispecies (1,2, and 
3) are from epidemiologically unrelated Americans 
5 and quasispecies 4 Is from a patient from 
Zimbabwe. 

Figure 6 shows polyacrylami.de heteroduplex 
mobility of intra- individual DNA heteroduplexes 
using complex DNA mixtures (quasispecies) as probe 
10 and target. Ml and M27 represent probes derived 
from the month 1 and 27 quasispecies. 

Figures 7A, 7B and 7C show polyacrylamide 
heteroduplex mobility of intra -individual DNA 
heteroduplexes using subclones derived from 
15 uncultured PBMC as probes and, quasispecies as 
target. 

Figure 8 shows a plot of the mobility shifts 
of heteroduplexes against the relative genetic 
distance between the two strands of the 
20 heteroduplexes. 

Figures 9A and 9B show the results of HMA- 
based and DNA sequence-based estimates of genetic 
relatedness for variants of the HIV-env gene. 

25 Detailed Description of the invention 

Definitions: 

"Homoduplex molecules" are typically composed 
of two complementary DNA strands, where the 
strands have at least about 98% sequence homology. 

30 "Heteroduplex molecules" are typically 

composed of two complementary DNA strands , where 
the strands have less than at least about 98% 
sequence homology. The functional definition of 
homoduplex and heteroduplex molecules, in the 

35 context of the present invention, is apparent from 
the results presented below. Typically, in a 
mixed population of homoduplex and heteroduplex 



14 

molecules, homoduplex molecules form a single band 
at a defined position in a polyacrylamide gel and 
heteroduplex molecules appear as slower migrating 
bands relative to the homoduplex. 

"Gaps" occur when in duplexes consisting of 
two complementary DNA strands, where th£ first 
strand of the DNA contains more nucleotides at an 
internal site than the second strand DNA molecule, 
and where these extra nucleotides are flanked by 
paired-complementary sequences . Gaps can occur in 
heteroduplexes. 

"Base-pair mismatches" typically refers to a 
single base-pair mismatch flanked by matched base- 
pairs. Base-pair mismatches also include a series 
of mismatched base-pairs flanked by matched base- 
pairs. Base-pair mismatches can occur in 
heteroduplexes. 

"Quasispecies" refers to the totality of 
members of a species whose genomes have sequence 
divergence, but where the members still fall 
within the same species: for example, outgrowth 
of "a number of sequence variants" (the 
quasispecies) of HIV- 1 from one major infecting 
variant. 

"Microorganism" in the context of the present 
invention includes, but is not limited to, the 
following groups: bacteria , viruses , fungi, 
parasites and mycoplasma. 

"Oncogene" includes genes that induce cancer 
or other uncontrolled proliferations of cells. 
Oncogenes can be mutated or activated (i) proto- 
oncogenes, or (ii) tumor suppressor genes (both of 
cellular origin) , associated with the development 
and/or proliferation of tumor cells. Oncogenies, 
or portions thereof, may also be of viral origin. 
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I. , . Observations on Hetei-n duplex Mobility in Non- 
Denaturina Pol van rvlamicte rc*» ip 

A - Altered Mobility of Heteroduplexes in 
Native Polvacryl.amide Gels. 

5 ■ Amplification reactions (Mullis; Mull is, et 
als) were performed on peripheral blood 
mononuclear cell (PBMC) DNA taken from an HIV 
seropositive asymptomatic man. The reactions were 
carried out using two sequential 25-35 primer 
10 extension cycles, where each cycle used a 
different set of primers (Example 1). The 
amplification products were a 1.8 kb (first round) 
fragment and then a 0.65 kb (second round) 
internal fragment which corresponded to an 
15 internal portion of the HIV-l envelope gene. When 
the amplification products were analyzed on an 
agarose gel, only the expected size band was 
observed (Example 2, Figure IB). However, when 
the amplified DNA was analyzed on a 5% 
20 polyacrylamide gel, additional, prominent bands of 
higher apparent molecular weight were observed 
(Example 2, Figure 1A) . 

Given the known variability of HIV-l within 
PBMC populations, experiments performed in support 
of the present invention suggested that the slower 
migrating bands may have been composed of 
heteroduplexes formed between divergent molecules 
during the last melt and reanneal (heat /cool) 
cycle. To explore this possibility, procedures 
30» expected to reduce or eliminate the formation. of 
heteroduplexes were performed. First, the 
heteroduplexes generated by PCR reactions were 
subjected to an additional single round of 
amplification (Example 2, Figure 1, lane 7). This 
35 procedure resulted in the loss of the slower 

migrating DNA bands. These amplification products 
were remelted and reannealed under conditions that 
prevented DNA synthesis, in this reaction 
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mixture, the series of slower migrating bands 
reappeared (Example 2, Figure 1, lane 8). 

Second, the concentrations of the initial 
PBMC DNA added to the amplification reactions was 
5 serially diluted (Example 2, Figure 1 , lanes 1- 
5) , This procedure also resulted in the loss of 
the slower migrating DNA bands. When the highest 
concentration of DNA was amplified and the 
resulting products denatured and renatured without 

10 amplification, the slower migrating bands were 
maintained (Figure 1A, lane 6) . 

In addition, if the additional bands were 
generated by heteroduplexes, then it seemed 
reasonable that such heteroduplexes/additional 

15 bands would be most efficiently formed at high 
product concentration, i.e.., at later rounds of 
amplification. As described above, amplification 
reactions were typically carried out with PBMC DNA 
as substrate and using two sequential 25-35 primer 

20 extension cycles was used to amplify a 1.8 kb and 
then a 0.65 kb internal fragment of the HIV- 1 
envelope gene. However, when only 25 cycles of 
amplification were used in the second round of 
amplification, the level of slower migrating DNA 

25 was less (Example 2, Figure 1, lane 9) than when 
more cycles of amplification were used. This is 
result is consistent with the idea that 
heteroduplex bands are most efficiently formed at 
high product concentration. 

30 All of the above-results are consistent with 

the conclusion that the additional bands in non- 
denaturing polyacrylamide gels were the result of 
heteroduplex formation during amplification 
reactions of HIV DNA from infected subjects. The 

35 results presented below confirm these 
observations. 
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B - The Effects of G aps and Substitutions on 
Heteroduplex Mobility, 

When two divergent sequences were mixed to 
form heteroduplex molecules, a single band was 
observed in agarose gels independent of the degree 
of sequence-relatedness between the two-strands of 
the heteroduplex (Figure 2 A) . However, in 
polyacrylamide gels, heteroduplexes formed from 
molecules with mismatched nucleotides, but without 
gaps, were detected when the degree of divergence 
reached above 1-2% and generally increased with 
the degree of mismatch (Figures 2B and 2D) . 
Mixtures of three different sequences yielded six 
extra bands (Figure 2C) . Thus each possible 
heteroduplex was formed* 

Further, heteroduplexes containing single 3 
bp internal gaps displayed mobility shifts (figure 
2E) . The effect of gapped sequences on 
heteroduplex mobility was further examined. A 
number of HIV-1 fragments from different source 
materials were amplified and fragments having 
divergent sequences were identified (Example 2, 
Figure 2E) . Based on sequence comparisons, three 
HIV-l fragments having internal deletions (i.e., 
deletions in HIV-1 sequences) of 3 and 9 base 
pairs were identified (3 and 9, Figure 2E) . 

Heteroduplexes were formed between the normal 
fragments, i.e., those not containing a deletion 
("insertion fragments", Figure 2E) , and the 
deletion fragments. The size of the gap present 
in the heteroduplex is shown across the top of the 
figure. These data show that centrally located 
gaps of 9 nucleotides resulted in heteroduplexes 
with slower migration than those with gaps of 3 
nucleotides. Further, these results indicate that 
the mobility of the heteroduplexes is affected by 
the sequence that is "looped out" of the 
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"insertion" sequence relative to the "deletion" 
sequence. 

The results presented above demonstrate that 
the sequence composition of the mismatches and 
5 gaps affects the mobility of the heteroduplexes in 
polyacrylamide gels. This observed effect may be 
adapted to probe for the presence or absence of 
specific nucleotides in pre-determined regions of 
interest, such as activating mutations in proto- 

10 oncogenes or particular resistance mutations in 
viral sequences. 

For example, since HMA could detect a 3 bp 
deletion in a 3.5 Kb DNA fragment it could also 
potentially be used to localize and utilize 

15 polymorphic sites for linkage mapping. Related 
micro-organisms could also be rapidly classified 
into sequence homology subgroups using HMA. HIV-1 
tissue culture and PCR contamination (problems of 
frequent occurrence) can be rapidly checked to 

20 ensure the identity and purity of the strain under 
study. 

C. The Relationship Between Heteroduplex 
Mobility Shift and DNA Sequence 
25 Distances. 

Phylogenetic relationships between genes 

encoding similar functions are usually described 

by the DNA distance expressed as the degree of 

similarity or mismatch between aligned sequences. 

30 These measurements ignore unpaired bases or gaps, 
because of the lack of a suitable method for 
assigning mutational distance arising from the 
latter types of mutations (Myers, G., et al. , 
1985a, 1985b). Having established that 

35 heteroduplex mobility shift occurs as a result of 
both base sequence divergence and gaps, the method 
of the present invention was used to evaluate the 
relationship between heteroduplex mobility and DNA 
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distance measurements based on DNA sequence 
• comparisons . 

The method of the present invention was 
employed to develop simple and rapid, assays of the 
5 relationships between members of DNA sequence 
families. In the following examples the 
heteroduplex mobility assays (HMA) are applied to 
the study of HIV-l envelope gene variability, but 
the method is also applicable to study of other 
10 viruses, microorganisms and gene families. 

DNA fragments (approximately 700 base pairs) 
encompassing the V3-V5 region of the HIV-l env 
gene from a set of 39 molecules of known sequence 
were amplified (Example 4) , These 39 molecules 
15 were as follows: thirty one molecules were 

obtained directly from the PBMC DNA of three North 
American subjects; six molecules from other 
independent HIV-l isolates from North America; and 
two HIV-l isolates from Zaire. Heteroduplexes 
20 were formed by pairwise combination of sequences 
from the same individual subject and between 
epidemiological ly unrelated North American and 
Zairian isolates. The relative mobility of the 
heteroduplexes compared to the known genetic 
25 divergence (Figure 3A) . Divergence was determined 
by standard methods of counting the number of 
mismatches between aligned sequences, discounting 
unpaired bases due to insertions and deletions 
(Kusumi, et al., 1992; Felsenstein, 1988). 

Three groupings were observed by both 
heteroduplex mobility and sequence analysis: (i) 
large mobility shifts and divergence were observed 
for heteroduplexes formed between Zairian/North 
American HIV-l strains; (ii) intermediate shifts 
and divergence were observed for heteroduplexes 
formed between independent North American 
isolates; and (iil) the smallest shifts and 
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divergence were observed for comparisons between 
the closely related sequences derived from the 
same subjects. 

Mobility shifts observed with some sequences 
5 derived from a long term infected asymptomatic 
individual (MA) were atypically large and 
partially overlapped those observed between 
independent US isolates (Figure 3A) . Further, the 
mobility shifts were generally lower in the HIV- 1 
10 pools from subjects with AIDS than from subject 
MA. When unpaired nucleotides were counted as 
mismatches the average genetic divergence 
increased more for the MA groups than for the AIDS 
patients. The relatively higher number of gaps 
15 distinguishing MA sequences is therefore likely to 
account for the greater mobility shifts observed. 
HIV-1 bjiv DNA derived from two MA PBMC samples 
collected 22 months apart were paired to each 
other for heteroduplex analysis. The mobility 
20 shifts of these heteroduplexes increased relative 
to heteroduplexes formed using the two DNA samples 
paired to DNA derived from the time point 
corresponding to when they were isolated. A 
similar increase in sequence divergence was also 
25 observed reflecting the evolution undergone by the 
virus pool during the sampling interval. These 
results suggest that the relative levels of HIV-1 
pool complexity (genetic divergence plus gaps) can 
be approximated by pairwise HMA of cloned 
30 sequences. 

Heteroduplex mobilities were plotted against 
genetic distance, as determined by percent 
sequence divergence (Figure 3A) . Heteroduplexes 
differing by only mismatches displayed a better 
35 correlation between mobility and divergence but 
with a lower slope (Figure 3A, circles) . The 
scatter of values in the intra-subject range may 
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be due to a large influence of . gaps on 
electrophoretic mobility. 

In other studies a more predictable 
relationship between mobility and DNA distance was 
5 observed between heteroduplexes formed between 
molecules twice as large, encompassing nearly the 
entire 1.3 kb extracellular envelope coding 
region. 

A standard curve for heteroduplex mobility 
10 can be generated using pair-wise combinations of 
DNA molecules having known genetic .distances, 
based on sequence comparisons between the 
molecules. For example, a region of a virus, such 
as the HIV env gene, can be selected, a number of 
15 variants are sequenced. Heteroduplexes are formed 
in pair-wise combinations between the variants. 
The genetic distance (i.e., percent sequence 
divergence based on mismatch) is then plotted 
against the heteroduplex mobility for each pair- 
20 wise combination. A sample of one such. plotted is 
presented as Figure 8. 

Heteroduplexes are then formed between pair- 
wise combinations of DNA samples, of the same 
region, by the method of the present invention. 
25 Heteroduplex mobilities are determined using HMA 
and the mobilities compared to the standard curve 
to establish the degree of relatedness between 
each pair-wise combination. This information ; can 
be used to generate phylogenetic trees such as is 
30 shown in Figure 9A. 

Figure 9A shows a phylogenetic tree based on 
HMA. Heteroduplex mobility data was recalculated 
into genetic distances using a HMA-DNA distance 
curve (Figure 8). The tree was generated using 
35 the "FITCH" program (Phylip 3.4 

shareware/ freeware) (Felsenstein, 1989). Figure 
9B shows a phylogenetic tree based on env gene DNA 
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sequences. Subtypes, identified by Myers (Myers, 
et al., 1992) are adjacent to boxes. Isolates 
present in Figure 9A that are not included in 
Figure 9B are indicated by an asterisk in Figure 
5 9A. 

Analyses such as this are useful for the 
determination of phylogenetic relationships 
between HIV-1 env genes or any other group of 
related sequences. 

10 

II. Applications 

The heteroduplex mobility assay (HMA) 
described herein provides a means for the rapid 
estimation of the degree of relatedness between 

15 members of gene families. HMA requires only a 

neutral polyacrylamide gel and yields results in 
hours. The method allows the determination of 
sequence relatedness without resorting to 
sequencing analysis (for example, Figure 9A) . 

20 Under the conditions described above, nucleotide 
gaps and mismatches can be observed in 
heteroduplexes using the method of the present 
invention. The preferred size range for these 
heteroduplexes is approximately 100 to 1500 base 

25 pairs in length, although larger-sized 
heteroduplexes can be used as well. 
Heteroduplexes formed from molecules with 
mismatched nucleotides, but without gaps, were 
detected when the degree of divergence exceeded 

30 approximately 1-2% and generally increased with 

the degree of mismatch. Typically, heteroduplexes 
containing an internal gap display mobility shifts 
which generally increased with the size of the 
gap. Further, 3 base gaps could be detected in 

35 heteroduplexes at least 3 kb in length. 

Accordingly, the method of the present invention 
provides the means to determine approximate levels 
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of DNA sequence diversity in a population of 
nucleic acid sequences both within and between 
individuals. The following examples using Human 
Immunodeficiency Virus sequence variations are 
illustrative of the method of the present 
invention. 

In addition to obtaining nucleic acid samples 
by amplification, other samples sources can be 
used as well. For example/ sequences of interest 
can be cloned (e.g.., in a lambda vector; Sambrook, 
et ai.) from two different sources. The sequences 
of interest are independently isolated away from 
vector sequences (e.g., by restriction 
endonuclease digestion and fragment purification) . 
These two samples can then be combined, denatured, 
renatured, and the resulting heteroduplexes 
analyzed as discussed below. 

A * DNA Sequence Div ergence Analysis. 

20 Estima tion of Degree of Sequence 

RelatednesB . 

Two previous studies have demonstrated 
outgrowth of a single genomic variant of env with 
a concomitant reduction in the diversity of HIV 
25 genes following the standard virus isolation 

method of in vitro co-culture with uninfected PBMC 
(Kusumi, et al., 1992). 

The method of the present invention was 
employed to compare the diversity of HIV-1 env 
30 genes found in PBMC versus those found after co- 
• culture. After co-culture, in 8/12 cases a single 
genomic variant, a descendant of an originally 
present variant, was observed. In 2/12 cases two 
genomic gap-phenotype variants were observed and 
35 in 2/12 cases three genomic gap-phenotype variants 
were observed, in each of the 12 independent 
cultures a marked reduction of diversity was 
observed. 
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When similar analysis was extended to 12 
other co-cuitures propagated from 2 to 8 weeks and 
a further 8 un-cultured in vivo quasispecies from 
different patients the pattern of low sequence 
5 complexity following co-culture and higher 
sequence complexity in vivo was maintained. 

In Figure 4 (Example 5), lane "Cult." 
presents HMA data for DNA samples amplified from a 
four week co-culture of HIV-infected PBMC. These 
10 results show a reduced sequence divergence in HIV- 
DNA obtained from co^cultured cells relative to 
the DNA derived from PMBC over time. Further, the 
comparison between two AIDS patients and an 
asymptomatic subject (Example 5) demonstrates that 
15 the overall level of HIV-1 env gene diversity in 

PBMC DNA (not co-cultured) samples was observed to 
be lower in subjects with low CD4 cell levels (PE 
and BU, Figure 4) than in those with higher CD4 
cell levels (MA time-points, Figure 4). 
20 The above results demonstrate that in vitro 

PBMC co-culture consistently results in the 
outgrowth of 1-3 distinguishable variants, 
confirming two previous studies based on extensive 
DNA sequence analysis (Meyerhans, et al., 1989; 
25 Kusumi, et al.; 1992), It is possible that the 
diverse viral genes observed in HIV-1 infected 
PBMC represent a record of prior infection in 
cells that survived an abortive infection with 
noninfectious virus. Alternatively, the culture 
30 milieu may impose selection pressure resulting in 
preferential replication of viral variants that 
are able to propagate quickly under the culture 
conditions. The heteroduplex mobility assay of 
the present invention may be used to distinguish 
35 between these possibilities: HMA can be used to 
effectively monitor the growth kinetics of 
specific variants in vitro. 
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The data presented above also demonstrate 
that eny gene complexity within the - blood of five 
subjects with AIDS or low CD4 cell numbers was 
lower than in six asymptomatic, seropositive 
5 subjects. Outgrowth of a, limited number of fast 
replicating variants in the presence of a 
declining immune system has been previously 
implicated in this apparent lower quasispecies 
complexity (Kuiken, et al., 1992; Tersmette, et 
10 al., 1989). 

2. Tracking se quence Variants . 
Experiments performed in support of the 
present invention have led to the development of 
an assay which rapidly estimates, DNA sequence 
variation without large scale DNA sequencing and 
which allows tracking and quantitation of variant 
genotypes in mixtures of such variants. Further, 
the method of the present invention allows such 
estimates with greater sensitivity than anything 
but very large scale DNA sequencing would allow. 
The heteroduplex mobility assay based on 
heteroduplex formation of PCR products is faster, 
simpler, and more informative than the currently 
available procedures (such as RNase A cleavage 
mismatch). HMA measures the retardation of 
electrophoretic mobility of DNA molecules in 
heteroduplex form relative to (matched) 
homoduplexes. 

The general method of the present invention 
has been described above. Following is a 
description of the application of the method of 
the present invention to the analysis of HIV-1 
envelope gene diversity, both within individuals, 
35 between individuals and between populations. In 
addition to analysis of viral genome diversity, 
such as described below for HIV-l, the method of 
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the present invention can be applied to the 
analysis of any number of microorganisms including 
bacteria, parasites, and other infectious agents. 
Exemplary microorganisms include, but are not 
5 limited to, the following: 

(i) Bacterial. Haemophilus — outer membrane 
proteins, Staphylococcus, Chlamydia — outer 
membrane proteins (Dean, et al.), Enterococcus , 
Mycobacterium (Mycobacterium tuberculosis) ; 

10 (ii) Viral. Feline Leukemia Virus (FeLV) , 

Simian Immunodeficiency Virus (SIV), Human 
Immunodeficiency Virus (HIV) , Hepatitis c Virus 
(HCV) ; Human papilloma virus (HPV) ; 

(iii) ^Fungi. Pneumococcus — Choline 

15 dependent Pneumococcal murein hydrolases; 18S rDNA 
sequences for human pathogenic fungi including 
Trichophyton, Histoplasma, blastomyces, 
coccidioides, Pneumocystis (Pneumocystis carinii) 
and Candida (Candida albicans) (Bowman, et al.); 

20 ( iv) Parasites . onchocerca ( 2 immerman , et 

al.Jr Babesia spp. (Ellis, et al.), Giardia spp. 
(Weiss, et al.), Leishmania spp. (Briones, et 
al.), Trypanosoma spp. (Breniere, et al.); and 
(v) Mycoplasma. Lyme disease, Mycoplasma 

25 pneumoniae (Kleemola, et al.)> using, for example, 
sequences derived from 16S RNA. 

Typically, probes for any target nucleic acid 
can be selected from a region of the 
microorganism's genomic material, such as rRNA 

30 (for example, as in Weisburg, et al.). In this 
way probes can be identified that will form 
homoduplexes to identify specific species. 
Formation of heteroduplexes indicates that the 
sequences that have diverged from the probe 

35 sequence. 

The method of the present invention can also 
be applied to the analysis of any nucleic acid 
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containing entity, including subcellular 
organelles such as chloroplasts and mitochondria. 

Further, the method of the present invention 
can also be used in screening methods for the 
evaluation of therapeutic treatments of any of the 
above microorganisms. The methods disclosed 
herein are useful for evaluating, in mixtures of 
nucleic acids (such as, nucleic acids obtained 
from tissue samples), the effect over time of a 
disease treatment, on DNA sequence variation of a 
nucleic acid target sequence associated with the 
disease. Therapeutic treatments typically are 
directed to the resolution, elimination, or relief 
of a disease state, as> for example, caused by a 
microorganism/ infectious agent. 

In one embodiment as applied to HIV 
infection, the heteroduplex mobility assay can be 
used to establish a base-line of infection in any 
selected patient before the onset of treatment. 
Diversity of the HIV virus can be determined as 
described below for different cell and tissue 
samples from the patient. Typically, blood and 
plasma samples are then serially collected from 
the subject throughout the therapeutic trials 

Changes during treatment in overall diversity 
of the HIV population can be monitored by the 
heteroduplex mobility assay. The turnover of 
specific variants can be monitored during the 
course of treatment use labeled tracer DNA (e.g., 
see below, "(b) Ouasisnecies Replacement within 
an Individual) . The emergence of new variants can 
also be detected by using specific HIV probes 
obtained from samples obtained from different time 
points in the course of the therapy. 
Representation of such variants in the population 
of HIV molecules can also be accomplished user 
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labeled tracer DNA (e*gr., Example 7, Tracking HIV- 
1 Sequence Variants V. 

In one embodiment, the method of the present 
invention is used to evaluate treatment of HIV 
5 with 3'-azidothymidine (AZT). The present method 
allows the evaluation of the effect of treatment 
on (i) quasispecies diversity, and (ii) the 
representation of particular variants within 
sample populations, in particular, variants pre- 

10 existing in the host that emerge after the onset 
of AZT therapy. Further, the method allows the 
determination of the time-frame in which any AZT- 
resistant mutants arise. during the course of 
treatment* After the identification of such 

15 mutants, the methods of the present invention can 
be used to assess the presence of such mutants in 
samples obtained before the mutants were observed 
(using, for example, the tracer method of Example 
7 ) . Further , the present methods can be used to 

20 follow the fate of such mutants during the course 
of further treatment; by, for example, 
dideoxyinosirie (DDI) . 

Tuberculosis (TB) infections are another 
example of how the method of the present invention 

25 can be used to monitor the effects of a disease 

treatment. The heterodup lex mobility assay can be 
used to monitor the presence and diviersity of 
strains of Mycobacterium tuberculosis growing 
within an individual. For example, a 383 bp 

30 segment of the gene encoding the 65 kDa 

mycobacterial surface antigen can be amplified 
(Ghossein, et al . , 1992) from samples obtained 
from a patient under treatment and analyzed by HMA 
and/or HTA. 

35 The assay can also be used to detect the 

specific loss or increase in abundance of TB 
variants during therapy. This is accomplished by 
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the tracer probe method of the present invention 
(Heteroduplex Tracking Assay (HTA) ; e.g., Example 
7) where the labeled probe is derived from 
standard TB strains or variants identified in the 
5 patient during the course of treatment. Such 
probes are used to track the representation of 
different Mycobacterium tuberculosis populations 
over time: before/ during and following therapy. 
The present methods can be used to help identify 
10 differential strain specificity to anti- 
tuberculosis drugs. 

Generally, the HTA methods of the present 
invention are used to monitor when variants come 
and go within the course of any infection and what 
15 the impact of any treatment has on the variant 

populations. Specific loci associated with drug 
resistance for a particular microorganisms can be 
used for tracking different populations of a 
microorganism using the methods of the present 
20 invention, where the variant loci are amenable to 
detection using HTA. 

The present assay can be used to evaluate 
diversity in cell culture systems and animal 
models as well as patients. 
25 Phylogenetic relationships can be established 

by the method of the present invention. 
Phylogenetic analysis can be carried out with 
almost any selected genomic sequence, such as, 
glycolytic enzymes (like phosphoglycerate kinase 
30 (Vohra, et al.)) or rRNA sequences. Phytogenic 
relationships between plants can be established, 
using, for example, sequences derived from plastid 
ribosomal RNA operons (Wolfe, et al.). 

Use of the method of the present invention to 
35 track sequence variants is described below using 
HIV-l as an exemplary microorganism . 
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(a) DNA Heteroduplex Analysis of Envelope 
Sequence Ouasispecies . 

Multiple (>50) HIV-1 envelope sequence 

variants were simultaneously amplified from 

5 peripheral blood mononuclear cell (PBMC) DNA by 

nested polymerase chain reactions (PCR) over the 

V3 to V5 regions of the env genes. The genetic 

relationship between different DNA quasispecies 

was determined by radioactively labeling one 

10 quasispecies DNA during PCR and reannealing it 
with an excess of unlabeled DNA from another 
quasispecies. Other labeling methods, in addition 
to radioactive labelling, can be employed: for 
example, labeling with biotin, digoxigenin, or 

15 chemi luminescent labels (Tropix, Inc. (Bedford, 
MA)). 

This procedure ensured that most of the probe 
tracer forms heteroduplexes with the unlabeled 
target driver DNA. The resulting heteroduplexes 

20 were then separated by polyacrylamide 

electrophoresis. When the probe and the target 
DNA quasispecies were identical, only fully 
complementary homoduplexes and moderately retarded 
radioactive heteroduplexes were seen reflecting 

25 the similarity between probe and target DNA, and 

the sequence heterogeneity within the quasispecies 
(Example 6, Figure 5) . 

When quasispecies from unrelated US sero- 
positives were similarly probed, no homoduplexes 

30 were seen and heteroduplexes migrated slower than 
did heteroduplexes from the same individual. When 
a quasispecies from a Zimbabwe sero-positive 
(Africa) was reannealed with US quasispecies, only 
very slowly migrating heteroduplexes were seen. 

35 The ability to probe one quasispecies with 

another allows a rapid determination of their 
degree of sequence similarity. Such information 
is valuable for epidemiological studies of 
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transmission risk factors by confirming or ruling 
out a source of infection (Ou, et al . , 1992). For 
determining the identity and purity of cultured 
isolates (e.g., HIV-1 strains) or PCR 
5 amplification products reactions: for example, if 
contamination is suspected (Palca, et al . ) . 

( b ) Quasispe cies Replacement within an 
Individual, 

10- Quasispegies change within an asymptomatic 

man was monitored by reannealing HIV-1 envelope 
DNA isolated at different time points spanning 27 
months. The progressive disappearance of labeled 
homoduplexes and the increasing mobility 

15 retardation of the labeled heteroduplexes as time 
increased between the probe and target 
quasispecies was observed. These results reflect 
the replacement of one set of sequences by another 
within two years (Figure 6, Example 6B) . The 

20 apparent sequence divergence between quasispecies 
separated by 27 months of in vivo evolution was 
still lower than that observed between 
epidemiological^ unrelated infected individuals 
from the US.' 

25 

(°) Tracking Sequence Variants in Comp lP.y 
DNA Mixtures. 

Sequential PMBC samples were obtained from an 
HIV-positive individual . Cloned variant sequences 

30 were used as probes to determine the prevalence of 
the cloned variants in the sequential PBMC 
samples. In these experiments the presence of 
radioactive homoduplexes indicated the presence of 
variant(s) highly related to the probe, i.e., the 

J5 cloned variants. Cloned variant sequences 

separated by 22 months of in vivo evolution were 
used to perform heteroduplex probing of six 
sequentially obtained PBMC samples collected over 
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27 months (Example 7) . Some variants were 
detected only in the quasispecies from which they 
were cloned (from independent PGR reactions) 
(Figure 7A, MA20 and Figure 7C, MA305) . Another 
5 cloned variant was found in PBMC samples obtained 
subsequent to the isolation of the cloned variant 
(Figure 7B, MA16) ; up to 22 months after the 
sample from which the cloned variant was isolated. 
A variant from the later time point was not 

10 detected in the preceding quasispecies (Figure 7C, 
MA305) . The patterns seen using molecular clones 
as probes, therefore, correspond to those seen 
using more complex probes (Figure 6, Example 6). 
However, use of the cloned probes provides a more 

15 detailed view of quasispecies change. 

The method of the present invention also 
provides a means to monitor the evo lilt ion of 
sequence complexity in cultured cells as well as 
in vivo. The method of the present invention has 

20 been used to track the presence of quasispecies of 
HIV-1 obtained from PMBC co-cultures. In addition 
to HIV-l tracking, the method of the present 
invention can be applied to variant tracking for 
other microorganisms. 

25 

(d) Evaluation of in vivo Quasispecies 
Complexity. 

The method of the present invention is also 
useful to analyze the in vivo evolution of 

30 quasispecies complexity for a given microorganism. 
PBMC DNA samples from four infected 
individuals prior to their sero-conversion were 
analyzed. From two samples, single molecules of 
HIV-1 DNA were amplified: 5 and 7 positive PCR 

35 end point sequences were obtained from these two 

samples. Pair-wise heteroduplex analysis of these 
PCR end points, using two different probes, showed 
no mobility shifts. This result indicates that 
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these pre-seroconversion quasispecies were made up 
of identical or very similar sequences. 

For two other quasispecies, a higher HIV-l 
DNA load allowed the simultaneous amplification of 
more than 100 envelope. sequences pre- 
seroconversion. In these samples the level of 
sequence complexity was directly estimated by 
polyacrylamide gel electrophoresis separation of 
its homoduplexes and heteroduplexes. This 
analysis showed these two pre-seroconversion 
quasispecies to be highly homogeneous. 

Similar sequence complexity analysis on 
sequential PBMC quasispecies was also performed. 
For two individuals (1058 and 527), the pre- 
seroconversion HIV-l DNA load in PBMC fell with 
the first signs of an antibody response to viral 
core antigens. Patient 537 HIV-l DNA developed 
slight heteroduplex mobility shifts 5.5 months 
following sero-conversion. Three months later, 
the quasispecies showed more substantial Variation 
concomitant with an increase in HIV-l DNA load. 

Patient 1058 first showed signs of low level 
envelope sequence complexity three months after 
the initial antigenemia peak. Further increase in 
quasispecies complexity was not apparent until 54 
months later, simultaneous with a noticeable 
increase in HIV-l DNA load and a detectable p24 
antigenemia. The result that HIV-l load (PBMC 
HIV-l DNA) fell with sero-conversion and then 
30 within months rose to intermediate levels is 

consistent with previous observations (Daar, et 
al.; Clark, et al.) . 

Highly homogeneous pre-seroconversion 
quasispecies DNA was used to probe subsequent 
35 quasispecies. Pre-seroconversion sequences were 
found in 537 PBMC at high levels for 14; 5 months 
post seroconversion and at reduced levels during 



20 



25 



WO 95/01453 



PCT/US94/07416 



34 

the next 32.5 months. At 54 months post- infection 
the quasispecies appeared free of the pre- 
seroconversion variant. 

Patient 1058> whose quasispecies appeared to 
5 remain homogeneous longer, retained sequences 
highly related to the major pre-seroconversion 
variant at high levels for the entire observation 
time of 57 months. A noticeable mobility shift 
beyond 21 months reflected the accumulation of 

10 nucleotide substitutions, which results in small 
mobility shifts, relative to the major infecting 
variant and the replacement of the quasispecies 
with such variants. Sequencing of PCR end points 
from each patient confirmed that (i) the pre- 

15 seroconversion quasispecies were highly related, 

and (ii) the later sequences had diverged from the 
pre-seroconversion sequences. 

3 . Geographical sequence variation . 

20 The method of the present invention also 

provides a means to track sequence variation 
between microorganisms isolated from different 
geographic locations by taking advantage of the 
gradual increase in mobility shifts observed with 

25 increasing sequence divergence between reannealed 
DNA strands. Analysis of this kind allows the 
rapid identification of sequence variants within 
species with, for example, viruses, such as, HIV 
and influenza. 

30 DNA fragments encoding the V1-V5 of the HIV-1 

envelope gene (env) from different geographic 
origins have been amplified and evaluated by the 
method of the present invention. Using 20 env 
fragments, whose DNA sequences were known, a plot 

35 was made of the heteroduplex mobility shift versus 
genetic distance. Prom the resulting curve, 
estimates were made of the genetic distance value 
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between unknown sequences based on the mobility 
shifts of their heteroduplexes. The genetic 
distance for all possible pairs of sequences from 
within the same country or continent was 
estimated. The sequences were as follows: ten 
sequences from Africa, eighteen from the 
US/Europe, six from Thailand and five from India. 
All of these sequences were derived from 
epidemiologically unrelated individuals. The 
mobilities of heteroduplexes made using HlV-i DNA 
found at the same time within uncultured PBMC from 
six individuals were also measured. 

The results of the mobility shift analyses 
were as follows. Thailand sequences fell into two 
.15 groups, m the first group the 4 independent 

isolates from Northern Thailand (Cheng-Mai) formed 
a very close cluster whose sequence diversity was 
comparable to that found within individuals . The 
two sequences from Southern Thailand (Bangkok) , 
20 while very divergent from those in Cheng-Mai, were 
themselves more related; 

The 5 independent Indian sequences from 
Bombay were significantly more related than 
independent sequences found in the US and Europe. 
25 African isolates exhibited high sequence 

divergence, although some Zambian/ Zimbabwe pair- 
wise combinations fell within the lower US/Europe 
diversity range. All pair-wise heteroduplex 
combinations tested using sequences within a 
30 subgroup resulted in electrophoretic migration 
typical of that subgroup (i.e., all 6 possible 
Cheng-Mai heteroduplex pairs showed a distinct 
grouping) , such that no heteroduplex pair behaved 
anomalously. 
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4. Phvlocrenetic Determinations bv 
Heteroduplex Analysis . 

To determine phylogenetic relationships 
between a number of isolates, heteroduplexes 
5 consisting of DNA strands from different sources, 
such as, different geographic, hospital or patient 
origins, can be analyzed by the method of the 
present invention. In the results described 
below, the first round primers were ED3 and ED14, 
10 and the second round primers were EDS and ED12 
(Example 1) . 

Such a phylogenetic analysis using the method 
of the present invention has been carried out as 
described above (Figures 8 and 9) . Subgroups of 
15 related sequences, transcending national and 

continental barriers, were apparent after mobility 
shift analysis. For some isolates, distant yet 
distinct relationships with particular subgroups 
were observed. Typically, the relatedness between 
20 groups was maintained between each pair tested 
(e.g., GD190, GD132, GD12 9 with the US/ Europe 
cluster; the Bombay, India with the major 
Zambia/ Zimbabwe clusters) . 

The subgrouping already apparent from a 
25 cursory examination of the mobility shifts was 
further refined using an algorithm to derive a 
phylogenetic tree based on heteroduplex mobility 
analysis (Felsenstein, et al., 1989). All 
examined US/European isolates clustered in a 
30 unique viral subgroup (Figure 9A) together with 
the two isolates from Bangkok (GD12 9, GD132) and 
one from Brazil (GD190) . A distinct subgroup was 
seen for 4 epidemiologically unrelated isolates 
from Somalia and Zambia and neighboring Zimbabwe 
35 (Zim 1, GD178, GD20, GD18-19 [sexual partners]). 
The sequence cluster from Bombay appeared to be 
related to this particular Eastern /Southern Africa 
cluster. 



Phylogenetic determinations of HIV-l isolates 
•from a number of geographical locations showed 
that the groupings obtained by heteroduplex 
mobility analysis correspond to those made with 
actual DNA sequences (Figure 9B) . Each of the 
"GD» HIV-l isolates studied here and provided 
under code were grouped in identical clusters by 
anchored PCR analysis in the gag region (samples 
were obtained from Dr. McCutchan; McCutchan, et 
al., 1992). It appears that HIV-l variation ' was 
greatest between isolates from the African 
continent, consistent with a longer residence time 
for the virus in Africa and larger resulting 
sequence diversification. 

A larger number of other Zambian strains were 
also reported to cluster using anchored PCR in the 
gag region (McCutchan, et al. , 1992). The 
Zambia/ Zimbabwe cluster were most related to the 
Bombay, India cluster. 

US/European sequences formed a large single 
related group with related sequences in Bangkok, 
Thailand and Brazil. A single HIV-l cluster in 
the US and Europe could reflect the descendence of 
most US/European isolates from one or a few highly 
related sequences in a manner similar to what is 
now seen in N. Thailand and Bombay. 

Independent isolates from Northern Thailand 
formed a very tight sequence cluster. The actual 
average nucleotide % substitution within four 
independent Thailand sequences for the V3 to V5 
region was 3;. 6% ranging from 2.6-4.5%, very close 
to what was found in four ihtra-patient 
quasispecies varying on average 2.3, 2.8, 3.3, and 
3.8% with a range of 0.15-8.41%. 

Three independent HIV-i sequences from 
Bombay, India also showed clustering with an 
average divergence of 6.1 % ranging from 4.8 to 
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7%. For contrast, 18 independent US sequence 
varied on average by 11% ranging from 7*6-13.8%. 

The ability to rapidly assign a microorganism 
(e.g., HIV) to groups of sequence homology should 
5 assist in determining the number of major 

subgroups in different locations. Since no tissue 
culture, subcloning, or sequencing is required for 
HMA, sequence variation can be rapidly determined 
by the method of the present invention. If 

10 successful vaccination is dependent upon or 
improved by a close match with the challenge 
strain the manufacture of vaccines related to 
local infectious microorganisms could increase 
vaccine efficacy. Vaccinated individuals, who 

15 nonetheless became infected, could be rapidly 

analyzed to determine the infecting microorganism 
subgroup. Such analysis should allow a more rapid 
appreciation of the effect of sequence variability 
on vaccine efficacy. 

20 Accordingly, the method of the present 

invention is useful for tracking quasispecies and 
sequence variations in samples containing viral 
nucleic acids. The method of the present 
invention can also be applied to tracking other 

25 pathogen-derived nucleic acids, for example, 
nucleic acids from bacteria, mycoplasma, 
protozoans, and parasites. Further, the results 
presented above demonstrate the usefulness of the 
method of the present invention for phylogenetic 

30 analysis and grouping of isolated microorganisms. 

This screening method can be applied to a 
number of microorganisms, as discussed above. The 
present method is useful in monitoring changing 
populations of gene sequences during the course of 
35 infection and disease. Further, the method is 
useful for tracking individual genomic variants 
and assessing overall levels of diversity, and for 



transmission of microorganism variants between 
individuals (e.g., in the case of HIV-1, maternal- 
fetal- transmission) and within populations (e.g., 
in the case of HIV-1 , in proposed WHO vaccine 
trial sites) . 

Homodupl ex Identification , 
Another embodiment of the present invention 
is the use of specific probes to identify variants 
based on the formation of homoduplex complexes. 
For example, sequences corresponding to a 
particular HIV variant can be cloned and 
amplified* These cloned sequences are then used 
as a probe against HIV molecules isolated from a 
number of test sources. Using the method of the 
present invention, if homoduplexes are formed in 
hybridization reactions between the probe and the 
test source HIV, then the test source HIV is shown 
to be similar to the cloned probe variant. If on 
the other hand heteroduplexes are formed between 
the probe and test sequences, then sequence 
divergence between the probe and test sequences is 
indicated. 

A battery of probe sequences can be 
established for major variants and then used in 
routine screening of test samples to identify the 
particular variant present in the test sample. 

As another example, probe sequences can be 
selected from the mycoplasma 16S rRNA for a number 
of species. Probe sequences for each species are 
generated by polymerase chain reaction using 
primer sequences conserved between mycoplasma 
species, but which flank regions of sequence 
divergence. DNA sequences corresponding to 
mycoplasma 16S rRNA sequences are available in the 
"GENBANK" data base for many mycoplasma species. 
A nucleic acid sample derived from a test source, 
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such as mycoplasma infected tissue culture, is 
then amplified with the same primers. The 
. resulting amplified test DNA is separately mixed 
with tracer amounts of labeled probe DNA (Example 
5 6) , denatured, renatured and resolved on 

polyacrylamide gels. Presence of homoduplexes 
indicates which species of mycoplasma is present 
in the infected source. 

This method is a self -confirming assay, in 
10 that, homoduplex formation should occur with only 
one probe sequence (representing one species) and 
heteroduplex formation with the other probes. 

C. Gapped Probes. 

15 Gapped probe molecules can be generated to 

regions of known sequence variation in a target 
gene. The mobility of heteroduplexes, formed 
between a target strand and a second strand 
containing an internal deletion relative to the 

20 target strand (i.e. < the gapped probe) , is 

affected by the sequence of looped out and the 
neighboring DNA sequences (Example 3, Figure 2E) . 
Accordingly, the mobility of such heteroduplexes 
can be used to evaluate the sequence present in 

25 the target molecule that bridges the gap. 

For example, the use of gapped probe 
molecules provides a way to identify the presence 
of mutations in a sample target strand relative to 
standard target strand that contains a wild-^type 

30 sequence. This method is particularly useful in 
the analysis of genes where the location of 
important mutations has been pr e-determined . 

D. Oncogene .Tracking. 

35 With respect to cancer, once a diagnosis has 

been made, and a region of DNA associated with the 
cancerous growth has been identified, the 
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heteroduplex tracking assay (HTA) of the present 
invention (e.g., Example' 7) can be used to 
evaluate the extent of infiltration of tumor cells 
within a tissue population, Exemplary potential 
target sequences are protooncogenes / for example, 
including but not limited to the following: c- 
myc, c-myb, c-fos, c-kit, ras, and BCR/ABL (e.g., 
Gazdar, et al., 1990; Wickstrom; Zalewski, et v al., 
1993; Calabretta, et al., 1992, 1993;), 
oncogenes/tumor suppressor genes (e.g., p53, 
Bayever, et al.). in tumor cells, deletions, 
insertions, rearrangements and divergent sequences 
in such genes or in the regions of DNA surrounding 
the coding sequences of such genes, all allow 
formation of heteroduplexes between amplified 
variant DNA and amplified DNA from normal cells. 

Specific probes can be designed to the 
variant oncogenic gene and the labeled probe can 
be used in HTA. The advantageous sensitivity of 
HTA is used to detect tumor cells within a 
population by heteroduplex tracking with 
sensitivities routinely as low as 1 in 100 cells. 
Thus, therapies that affect his prevalence can be 
ascertained through the use of HTA. 

The following examples illustrate, but in no 
way are intended to limit the present invention. 

Materials and Methods 
E. coli DNA polymerase I (Klenow fragment) 
was obtained from Boehringer Mannheim Biochemicals 
(Indianapolis, IN) . T4 DNA ligase and T4 DNA 
polymerase were obtained from New England Biolabs 
(Beverly, MA) ; Nitrocellulose filters were 
obtained from Schleicher and Schuell (Keene, NH) . 
Restriction enzymes were purchased from commercial 
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vendors and used as per the manufacturer's 
instructions. 

Synthetic oligonucleotide linkers and primers 
were prepared using commercially available 
5 automated oligonucleotide synthesizers* 

Alternatively, custom designed synthetic oligo- 
nucleotides may be purchased, for example, from 
Synthetic Genetics (San Diego, CA) . DNA labeling 
kits can be obtained from Bdehringer-Mannheim 
10 Biochemicals (BMB, Indianapolis, IN) and Bethesda 
Research Laboratories (Gaithersburg, MD) . 
Polymerase chain reaction reagents and equipment 
are available from Perkin Elmer Corporation 
(Norwalk, CT) . 

15 Routine molecular biology manipulations were 

carried out by standard procedures as taught, for 
example, by Ausubel, et al., Mianiatis, et al., or 
Sambrook , et al. 

20 EXAMPLE 1 

Heteroduolex Formation and Analysis 

A. Cellular DNA Extraction and Polymerase 
Chain Reactions * 

Peripheral blood mononuclear cells (PBMC) 

25 isolated by "FICOLL/HYPAQUE" density gradient 

centrifugation were washed twice with phosphate 

buffered saline (PBS; Gibco-BRL, Gaithersburg MD) . 

DNA was extracted from the PBMC using the 

"ISOQUICK DNA" isolation kit (MicroProbe Corp., 

30 Garden Grove, CA) . 

Polymerase chain reactions (PCR) (Mullis; 

Mullis, et al.) were carried out in two rounds 

using nested primers. Typically, 2 ul of the 

first round reaction product was added to a second 

35 round of PCR with internally annealing primers. 

First round primers were typically selected from 

one of the following groups: ED3 (SEQ ID NO:l), 

corresponding to positions 5956*5985 on the HIV-1- 
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HXB2 genome (Ratner , et al . ) and ED12 (SEQ ID 
NO: 2), corresponding to the complement of 
positions 7810-7781 on the HIV-1-HXB2 genome; or 
ED3 and ED14 (SEQ ID NO: 7) , corresponding to 
5 positions 7936-7966 on the HIV-1-HXB2 genome. 

One set of second round primers was ES7 (SEQ 
ID NO: 3), corresponding to the M13 universal 
primer followed by positions 7000-7020 of HIV-1- 
HXB2) and ESS (SEQ ID NO: 4) , corresponding to the 

10 complement of the M13 reverse sequencing primer 
followed by positions 7667-7647 of HIV-1-HXB2) . 
The second round primers give rise to an 
amplification product of 704 bp of which 
approximately 627 bp are template dependent: the 

15 actual size of the amplification products depends 
on the size of deletions and insertions within the 
target molecule relative to the HIV-1-HXB2 
template. 

A second set of second round amplification 

20 primers contain the sequences ED5 (SEQ ID NO: 8), 
positions 6562-6588 of HIV-1-HXB2, and ED12 (SEQ 
ID NO:2), complement of positions 7792-7822 of 
HIV-1-HXB2. For the geographic analysis, first 
round primers contained ED3 and ED14 and second 

25 round primers contained ED5 and ED12 . 

A third set of second round amplification 
primers were ESS (SEQ ID NO: 5), positions 7521- 
7540 of HIV-1-HXB2, and ED4 (SEQ ID NO:6), 
complement of positions 7741-7712 of HIV-1-HXB2. 

30 Use of primers ESS and ED4 yielded an -220 bp 
amplification fragment. 

Each PCR reaction employed variable amounts 
of template DNA (up to 1 fig) , 1.8 mM MgCl2, 20 
pmole of each primer in 50 mM KC1, 10 mM Tris-HCl 

35 pH8.3, 200 uM of each dNTP, 2.5 units of Taq DNA 
polymerase (Perkin Elmer-Cetus, Emeryville, CA) , 
and 10% glycerol in a final volume of 50ul. PCR 
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reactions were carried out, essentially as per the 
manufacturer's suggestion, in a perkin Elmer 
Thermocycler Model 124 for 25-35 cycles using 1 
sec ramp times between steps of 94 °C for 1 minute, 
5 57°C for 45 sec, and 72°C for 1 minute. 

For heteroduplex resolution a 5 minute 72 °C 
extension step was linked to the last cycle, 
Heteroduplex resolution was carried out by 
transferring 10 ill of the second round PCR 
10 reaction to 90 jil of fresh standard PGR reaction 
mix (with 100 pmole of each primer) followed by 
one denaturation and 5 minute extension cycle at 
72°C. 

15 B. End point dilution determination of 

viral DNA load . 

Titration of HIV-1 DNA in PBMC DNA was 

performed essentially as described by Simmonds, et 

al., (1990). Briefly, HIV-1 DNA was titrated 

20 using duplicate serial 5 fold dilutions of input 
infected PBMC DNA and maintaining a constant 1 ug 
of human genomic DNA. 

The lowest concentration of infected cell DNA 
to consistently yield a positive HIV-1 PCR signal 

25 was used to estimate the proviral DNA load. When 
a low proviral DNA load precluded simultaneous 
amplification of at least 20 template molecules, 
products from multiple reactions, were pooled. As 
controls to demonstrate single molecule template 

30 sensitivity, pNL4-3 (HIV-1 Lai) DNA (Adachi, et 
al.) and ACH-2 cell DNA (containing a single 
defective HIV-1 genome) (Folks, et al.) were used 
as infected DNA sources and diluted in the 
presence of a total of 1 of human genomic DNA. 

35 Single molecule HIV DNA templates from 

subject PBMC were also derived by endpoint 
dilution and used to generate probe sequences for 
tracking their representation in PBMC DNA over 



time, Endpoint ; ^reactions were subjected to the 
heteroduplex mobility assay of the present 
invention to verify template homogeneity. 

PCR and cloning artifacts were minimized by 
sequencing HIV-1 variants directly from PCR 
fragments derived by nested PCR from single 
molecules of HIV-l (Simmonds, et al . , (1990)). 
This procedure also provides variant frequencies 
representative of the prevalence of the HIVr-i 
molecules within the quasispecies; as may not 
happen following subcloning, when a single 
provirus can give rise to two or more subclones. 

DNA sequencing of isolated amplification 
products or. clones of HIV DNA was carried out by 
standard dideoxy sequencing reactions (U.S. 
Biochemical, Cleveland OH) . 

c - Heteroduplex formation and analysis . 
Typically, 4.5 /xl of PCR product from 2 
separate reactions were combined and l-l . 5 /il of 
10X annealing buffer added (1M NaCl^ lOOmM Tris 
HC1 pH 7 .8, 20 mM EDTA) . The mixtures were heated 
to 94°C for 2-5 minutes and rapidly cooled to 22°C 
in the thermocycler (1 second ramp). For analysis 
of HIV samples with high sequence divergence 
(e.g., samples from different geographical 
locations) the mixtures were heated to 94 6 C for 2- 
5 minutes and rapidly by placing the samples in an 
ice bath. ...... : • 

Heteroduplexes were then separated on 5% 
polyacrylamide gels (30:0.8 acrylamiderBis) at 
250V for 3 hours in IX TBE (0.088M Tris-borate, 
0.089M boric acid, 0.002M EDTA) or 2.5% agarose 
gels at 100V for 1.5 hrs. in TAE buffer (0.04M 
Tris-acetate, 0.001M EDTA). The gels were stained 
with ethidium bromide and photographed. The 
temperature at which the gels are maintained 



during electrophoresis can affect the mobility of 
the heteroduplexes. 

EXAMPLE 2 

Reduced Mobility of 650 bp DNA H eteroduplexes 
in Native Polvacrylamide Gels 

Nested PCR using two sequential 25-35 primer 

extension cycles was used to amplify a 1.8 kb and 

then a 0.65 kb internal fragment of the HIV-1 

envelope gene directly from peripheral blood 

mononuclear cell (PBMC) DNA taken from an HIV 

seropositive asymptomatic man (Example 1). When 

the sample was analyzed on a 2.5% agarose gel, 

only the expected size band was observed (Figure 

IB, lanes 1-5) . However, when the DNA was 

analyzed on a 5% polyacrylamide gel (Example 1) , 

additional, prominent bands of higher apparent 

molecular weight were observed (Figure 1A, lanes 

1-5) . 

To determine the nature of these additional 
bands, the following analyses were carried out. 
First, a fraction of the PCR reaction was 
subjected to a single additional round of PCR 
using fresh polymerase and an excess of primers 
(Figure 1, "Resolve", lane 7). Second, the 
concentration of PBMC DNA was serially reduced in 
the PCR reactions (Figure 1, lanes 1-5), 

Each of these manipulations resulted in the 
loss of the slower migrating DNA bands (Figure 
1A) . This result supports the conclusion that the 
additional bands seen in the initial PCR reactions 
were heteroduplexes formed between divergent 
molecules during the last melt and reanneal 
(heat /cool) cycle of the PCR reactions. 

Figure 1A, lane 6 shows the result of melting 
and reannealing an aliquot of the sample that was 
fractionated in lane 5. This melting and 
annealing was not accompanied by an additional 
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round of amplification. As expected, the number 
•of bands in each of lanes 5 and 6 (Figure 1A) are 
essentially the same. 

Furthermore, the DNA which had been subjected 
to a single additional round of PCR using fresh 
polymerase and an excess of primers- (sample used 
in Figure l, lane 7) was remelted and reanriealed 
in the presence of EDTA to prevent Tag polymerase 
activity. When this sample was resolved on a 5% 
polyacrylamide gel, the series of slower migrating 
bands reappeared (Figure 1A, Heat/Cool, lane 8). 

In addition, nested PCR was carried out using 
only one round of 25 primer extension cycles to 
amplify a 0.65 kb internal fragment of the HIV-l 
envelope gene directly from the PBMC DNA-described 
above, in these reactions, the products 
demonstrate that the level of slower migrating 
DNA, i.e., the additional bands, was reduced 
(Figure 1, lane 9) relative to the 35 cycle 
20 amplification described above. 

Each of the results presented above are 
consistent with heteroduplex formation during 
later stages of amplification when HIV DNA from an 
infected subject is used as template DNA. 



15 



25 



EXAMPLE 3 

Heterodunlevft s Containing Sequence Gap s 

or Base-P air Mismatches Demonstrate 
a Mobility Shift in Polyacrylamide Gels 

30 A « Mismatch and Gap s. 

Heteroduplexes were formed by melting and 
reannealing mixtures of DNA fragments essentially 
as described above. Samples of HIV-l were 
obtained from a variety of sources. Briefly, 

35 molecular clones of HIV-l genes in plasmids were 
obtained (e.g., Kusumi, et al.). Sequences were 
either obtained from published sequences or by 
employing standard sequencing techniques with 
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using isolated HIV sequences. Heteroduplexes Were 
formed as described in Example l. The migration 
of the heteroduplexes and homoduplexes were 
compared in polyacrylamide and agarose gels as 
5 described in Examples 1C and 2. 

Figures 2B and 2C show photographs of 
ethidium bromide stained 5% polyacrylamide gels 
containing, in each lane, heteroduplexes formed 
from a mixture of two PCR products, each PCR 

10 amplified product has been sequenced. The level 
of sequence diversity between the two products is 
indicated in the figure: sesquence percent 
mismatch ranges from 0.16 to l.ll. Also, the 
presence (+) or absence (-) of gaps in the 

15 heteroduplex sequences are indicated. 

In Figure 2D results are shown for 
heteroduplexes having the levels of sequence 
diversity indicated in the figure: percent 
mismatch from 1.3 to 4.9. 

20 When two divergent sequences were mixed, a 

single band was observed in agarose gels 
regardless of the sequence relationships between 
the sequences (Figure 2A) . In polyacrylamide gels 
the same mixtures resulted in nearly comigrating 

25 homoduplex bands (bottom band in each lane, 

Figures 2B and 2C) plus two additional slower 
migrating bands (Figures 2B and 2C) . 

Mixing and annealing three different-sequence 
amplification products yielded six bands (Figure 

30 2D, lane 1) in addition to the homoduplex band 

(which corresponds to the fastest migrating band 
in Figure 2D, lane 1). Thus each possible 
heteroduplex is formed. Figure 2D, lane 2 
contains molecular weight markers derived from a 

35 Haelll digest of 0X174. 

The results presented above suggest that the 
composition of the mismatches and gaps affects 
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mobility. The effect of different sized gaps on 
heteroduplex mobility was next examined. 

B - The Effects of se quence Variations 1 
Within Insertion SeauennBs. 

A number of HIV-1 fragments from different 
source materials were amplified (second round 
primers ES7 and ES8) and sequenced as described 
above (Example 1). The amplification products 
were normally 704 bp of which approximately 627 bp 
are template dependent. A number of 704 bp 
fragments having divergent sequences were 
identified ("Insertion" fragments, A-H, Figure 
2E). The duplicate lanes correspond to fragments 
15 having the same nucleotide sequence that were 

derived from different sources. Further, based on 
sequence comparisons, three HIV-l fragments having 
internal deletions (i.e., deletions in HIV-1 
sequences) of 3 and 9 base pairs were identified 
20 (3 and 9, Figure 2E) , which were designated MA21 

(3 base pairs (bp)), MA311 (9 bp), and MA 6 (9 bp). 

Heteroduplexes were formed between the 
insertion fragments and the deletion fragments. 
The heteroduplexes were electrophoretically 
25 separated as described in Example l. The gels 
were stained with ethidium bromide and 
photographed. The results are shown in Figure 2E. 
In the figure, the bottom band in each lane 
corresponds to homoduplexes . The lanes marked "M" 
30 are molecular weight standards. 

The size of the gap present in the 
heteroduplex is shown across the top of the 
figure. These data show that centrally located 
gaps of 9 nucleotides resulted in heteroduplexes 
35 with slower migration than those with gaps of 3 

nucleotides. Further, these results indicate that 
the mobility of the heteroduplexes is affected by 
the sequence that is "looped out" of the 
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"insertion" sequence relative to the "deletion 11 
sequence. For example, heteroduplexes formed with 
deletion fragment 12 have a range of mobilities 
wher;e the mobilities are dependent on the 
5 sequences present in the looped out insertion of 
the second fragment A, B, C or D. 

EXAMPLE 4 

Heteroduplex Mobility Shift and 
10 DNA Sequence Distances 

A number of HIV-1 fragments, corresponding to 

HIV-1 env sequences (V3-V5) from within the same 

and between different seropositive individuals, 

were amplified (second round primers ES7 and ES8) 

15 and sequenced as described above (Example 1) . The 
amplification products were normally 704 bp of 
which approximately 627 bp are template dependent. 

Relative mobilities of HIV-1 DNA intra- and 
inter-subject heteroduplexes, formed using 

20 amplified DNA from the following sources, were 

evaluated by the method of the present invention: 
MA89-91 is a comparison of sequences derived from 
subject MA 22 months apart in 1989 and 1991; US/US 
is a comparison of sequences from 

25 epidemiologically distinct viruses from subjects 

within the United States of America; US/AFR is a 
comparison of sequences from 6 US [pNL4-3, SF2, 
SF162, MA5 (Myers, et al., 1992); BU01, PE01] and 
2 Zairian (Africa) subjects [NDK> MAL (Myers, et 

30 al., 1992)]. Included in the US/AFR group is the 
NDK/MAL comparison which displayed the fastest 
mobility and least sequence divergence in the 
group. 

For each comparison of heteroduplexes formed 
35 using amplification products from two source DNAs, 
heteroduplex mobilities on non-denaturing 
polyacrylamide gels were calculated as the average 
distance of migration of the two heteroduplex 
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bands divided by the distance of migration of the 
homoduplex bands. The sequence divergence between 
the nucleic acid sequences from each HIV-fragment 
source were determined by comparison of the 
nucleic acid sequences using the program DOTS 
(Kusumi, et al., 1992). The program counts the 
number of mismatched bases between aligned 
sequences, discounting gaps introduced to maintain 
alignment. 

In Figure 3A the relative mobility values are 
plotted against the percent divergence for four 
sets of sequence comparisons: +, Intrasubject, 
from within, subjects MA, PE and BU; • , 
Intrasubject, including only comparisons for which 
15 no unpaired segments or gaps appear within 

heteroduplexes; n, us/US inter-subject; O, US/AFR 
inter-subject. 

Figure 3B shows representative heteroduplex 
mobility data using the above DNA sources and a 5% 
20 polyacrylamide gel with representative 

heteroduplexes formed from the intrasubject, 
intersubject and US/AFR groups, respectively. 

The above data demonstrate the correlation 
between the degree of sequence relatedness, based 
on direct nucleic acid sequence comparisons, and 
heteroduplex/homoduplex mobility on non-denaturing 
polyacrylamide gels: generally,, the relative 
mobility shift increases with increasing percent 
divergence within patient samples. 



25 



30 



35 



EXAMPLE 5 

Consistent Reduc tion of Pool Diversity 
Upon PBMC-Co-Cultiiyo 

The heteroduplex mobility assay described 
above was used to compare the diversity of HIV-l 
env genes found in isolated PBMG versus those 
found after co-culture of the PMBC (Kusumi, et 
al., 1992). 



52 

DNA was isolated from PMBC as described in 
Example 1. Titration of the HIV-1 DNA in the PBMC 
DNA was performed as described in Example 1. Each 
PBMC DNA sample contained at least 20 molecules of 
HIV-1 DNA, determined by end point PCR. In each 
sample the HIV-1 molecules were simultaneously PCR 
amplified as described above using ES7 and ES8 as 
second round primers in nested PCR reactions. 

Figure 4 shows the results of heteroduplex 
mobility analysis of viral quasispecies obtained 
from PBMC and after virus isolation in co-culture. 
In the figure, lane numbers refer to the month of 
collection of PBMC from an asymptomatic subject 
(MA) beginning approximately 5 years after 
infection (Kusumi, et al., 1992)* Duplicate 
nested PCR reactions are shown for the 1 and 7 
month time points. "Cult." refers to DNA samples 
amplified from a four week co-culture of HIV- 
infected PBMC (time point, 1 month) with 
uninfected PBMC. Lane W M M contains molecular 
weight markers. "pe m and "BU M correspond to 
source PMBC DNA obtained from two AIDS patients. 

The results presented in Figure 4 indicate a 
reduced sequence divergence in HIV-DNA obtained 
from co-cultured cells relative to the DNA derived 
from PMBC over time. Further, the PMBC DNA 
samples from the two AIDs patients, PE and BU, 
also demonstrate a lower sequence divergence 
relative to the amplified samples from the 
asymptomatic patient. 

EXAMPLE 6 

Radioactive Probes and HMA 

A. Relationships between Quasispecies from 
Epidemiological^ unlinked Individuals. 

Typically, multiple (>50) HIV-1 envelope 

sequence variants were simultaneously amplified 

(Example 1) from peripheral blood mononuclear cell 
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(PBMC) DNA by nested polymerase chain reactions 
(PGR) over the V3 to V5 regions of the env genes 
(Example 1, using ES7 and ES8 as second round 
primers in nested PCR reactions) . 

The genetic relationship between DNA 
molecules amplified from different source DNAs was 
determined by radioactively labeling one group of 
DNA molecules during PCR and reannealing those 
labeled molecules with an 100-fold excess of 
unlabeled DNA from another source material. 

With this procedure, most or all of the probe 
"tracer" (i.e., labelled amplified DNA molecules) 
formed heteroduplexes with the unlabeled target 
"driver" DNA. The resulting heteroduplexes were 
then separated by polyacrylamide electrophoresis 
(Example 1) . The gels were then exposed to X-ray 
film (with or without intensifying screens) and 
the resulting autoradiograms analyzed. 

Figure 5 shows the results of one such 
"tracer"/ "driver" analysis. PMBC DNA was isolated 
from four sources, three US HIV-positive samples 
and one African HIV-positive sample (from 
Zimbabwe) . Nested DNA amplification reactions 
were carried out as described in Example 1 using 
primers ES7 and ES8 . 

The probes were labeled PCR products from 
each PMBC sample. The products were radiolabeled 
by addition of 10 pci of a 32 P-dCTP and 30 of 
each dNTP to the second round of a nested PCR 
reaction. 

Unlabeled PCR products from the subject PBMC 
samples were mixed, heated and reannealed in 100 
fold excess with each of the radiolabeled probes. 
The resulting reannealed products were 
electrophoretically separated on a 5% 
polyacrylamide gel. An autoradiograph of four 
such analyses is shown in Figure 5. The asterisk 
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(*) to the left of the panel denotes the position 
of single stranded DNA, which is of variable 
intensity and pattern from experiment to 
experiment* 

5 When the probe and the target DNA were the 

same (Figure 5, lanes 1-1 ,, 2-2, 3-3, and 4-4), 
only fully complementary homoduplexes and 
moderately retarded radioactive heteroduplexes 
(representing the members of this quasispecies) 
10 were seen. This result reflects the similarity 
between probe and target DNA, and the limited 
sequence heterogeneity within the quasispecies* 
When quasispecies from unrelated US sero- 
positives were similarly probed, no homoduplexes 
15 were seen and heteroduplexes migrated slower than 
did heteroduplexes from the same individual 
(Figure 5, panels 1-3, lanes 1-3). When a 
quasispecies from a Zimbabwe sero-positive 
(Africa) was reannealed with the US quasispecies, 
20 only very slowly migrating heteroduplexes were 

seen (Figure 5, panels 1-3, lane 4> and panel 4, 
lanes 1-3) . 

B. Intra-Individual HMA for HIV-1 
Quasispecies * 

Quasispecies change within an asymptomatic 
HIV-positive man (MA) was monitored as follows. 
HIV-l envelope DNA was isolated by nested PCR 
amplification, using ES7 and ES8 as second round 
primers, from PMBC DNA. The PBMC DNA was isolated 
from the subject at different time points spanning 
27 months. "Tracer" /"driver" analysis was carried 
out as described in Example 5 using amplified DNA 
from each PMBC time point. 

In this example, the source of the PMBC DNA 
was asymptomatic patient MA, described above in 
Example 5. 



25 
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Probes were derived- from PGR products from 
the month 1 (Ml, Figure 5, left panel) or month 27 
(M27, Figure 5, right panel) j Month 1 and 27 
products were radiolabeled by addition of 10 [id 
5 of a J2 P-dCTP and 30 juM of each dNTP to the second 
round of the nested PCR reaction . 

Polymerase chain reaction products from the 
asymptomatic subject MA PBMC samples, described 
above in Example 5, were mixed, heated and 
10 reannealed in 100 fold excess with two 
radiolabeled probes . 

The resulting reannealed duplex products were 
electrophoretically separated on a 5% 
polyacrylamide gel. An autoradiograph of the gel 
15 is shown in Figure 6. The asterisk (*) to the 
left of the panel denotes the position of single 
stranded DNA. 

The labeled probe DNA (tracer) is shown at 
the top of the figure: Ml corresponds to the 
20 quasispecies amplified from the first month PMBC 
and M27 corresponds to the quasispecies amplified 
from the 27 month PMBC. The driver DNA is shown 
in the second line of the figure: the number 
represents the month the sample was obtained. 
55 The progressive disappearance of labeled 

homoduplexes (Ml-tracer) and the increasing 
mobility retardation of the labeled heteroduplexes 
(Ml-tracer) as time increased between the probe 
and target quasispecies was observed. The M27- 
to tracer demonstrates the replacement of one set of 
sequences (M27-tracer :Ml-driver) by another (M27- 
tracer:M2 7 -driver) within two years. 

EXAMPLE 7 

5 Tracking HIV-T sequence Variants 

PMBC DNA was isolated from an HIV-positive 
individual over a twenty seven month period. 
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Variant sequences were cloned from the 1 month 
PMBC sample (clones MA20 and MA16) and from the 22 
month sample (clone MA305) : the source DNAs for 
these clones were the unamplified PMBC DNA 
5 samples. HIV sequences (for example, env 

sequences) were cloned using standard vectors and 
techniques (Kusumi, et al.> 1992; Ausubel, et al.; 
Sambrook, et al.). The cloned sequences were 
labeled with radioactive moieties by end-labeling, 
10 random-priming or nick-translation (kits for each 
method available from Gibco/BRL, Gaithersburg, 
MD). 

The HIV eiiv sequences were amplified as 
described in Example 1 using primers ES7 and ES8 

15 and the sequentially obtained PBMC DNAs as 

template (as in Example 6) • The cloned DNAs were 
used as probes to determine their prevalence in 
sequential PBMC samples as described in Example 6. 
The results of this analysis are shown in Figures 

20 7A, 7B, and 7C. 

In Figure 7, the numbers at the top of the 
lanes correspond to the month that the sequential 
PMBC sample was isolated. The clone names are 
indicated at the top of each panel of the figure 

25 (Figure 7A is MA20, Figure 7B is MA16, and Figure 
7C is MA305) . Ml indicates that the probe was 
obtained from the 1 month PMBC DNA sample. M22 
indicates that the probe was obtained from the 22 
month PMBC DNA sample, "s.s." is the location of 

30 the labeled single-strand DNA and "H" is the 
location of the homoduplex. 

One clone (MA20) obtained from the 1 month 
PMBC DNA detected the presence of its 
corresponding variant only in the quasispecies 
35 (i.e., the family of variant HIV sequences present 
in the PMBC sample) from which it was cloned ("H" , 
lane 1, Figure 7A) . Another clone, MA305, which 
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represents a different quasispecies, was present 
in the amplified PMBC DNA from which it was 
obtained (Figure 7C, lane 22) and also in two 
subsequent samples (Figure 7C, lanes 23 and 27) . 
5 The third clone, MA16, was found in the 

amplified PMBC DNA from which it was obtained 
(Figure 7B, lane l) and also in subsequent PBMC 
samples up to 22 months later (Figure 7B, lanes 6, 
7 and 22) . 

10 In Figures 10B and IOC, the lane marked »C M 

is a control reaction showing the results of 
mixing the tracer clone DNA with the driver DNA 
corresponding to the PMBC sample from which the 
clone was obtained. 

15 

While the invention has been described with 
reference to specific methods and embodiments, it 
will be appreciated that various modifications and 
changes may be made without departing from the 
20 invention. 
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SEQUENCE LISTING 

(1) GENERAL INFORMATION: 

5 (i) APPLICANT: 

(A) NAME: The Board of Trustees of the Leland Stanford 
Junior University 

(C) CITY: Stanford 

(D) STATE: CA 
10 {£) COUNTRY: USA 

(P) POSTAL CODE: 94305 

(i) APPLICANT: 

(A) NAME: The Government of the United States of America 
15 as represented by The Secretary of the Department 

of Health and Human Services and His Successors 

(C) CITY: Washington 

(D) STATE: DC 

(E) COUNTRY: USA 

20 (P) POSTAL CODE: 20231 

(ii) TITLE OP INVENTION: A Heter ©duplex Mobility Assay for the 
Analysis of Nucleic Acid Sequence Diversity 

25 (iii) NUMBER OF SEQUENCES : 8 

(iv) CORRESPONDENCE ADDRESS: 

(A) ADDRESSEE: Dehlinger and Associates 
. (B) STREET: 350 Cambridge Avenue, Suite 250 
30 (C) CITY: Palo Alto 

(D) STATE: CA 

(E) COUNTRY: USA 

(F) ZIP: 94306 

35 (v) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS /MS-DOS 

(D) SOFTWARE: Patent In Release #1*0, Version #1.25 

40 

(vi) CURRENT APPLICATION DATA: 

(A) APPLICATION NUMBER: 

(B) FILING DATE: 

(C) CLASSIFICATION: 

45 
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(vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: US 08/241,373 

(B) FILING DATE: ll-MAY-1994 

5 (vii) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: US 08/087,010 

(B) FILING DATE: l-JUL-1993 

(viii) ATTORNEY /AGENT INFORMATION: 
10 (A) NAME: Fabian, Gary R. 

(B) REGISTRATION NUMBER: 33,875 

(C) REFERENCE /DOCKET NUMBER: 8600-0130.41 

(ix) TELECOMMUNICATION INFORMATION: 
15 (A) TELEPHONE: (.415) - 324-0880- 

(B) TELEFAX: (415) 324-0960 

(2) INFORMATION FOR SEQ ID NO:l: 

20 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDED NESS : single 

(D) TOPOLOGY: linear 

25 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

30 (vi) ORIGINAL SOURCE: 

(C) INDIVIDUAL ISOLATE: PRIMER BD3 , HIV-1-HXB2 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:l: 

35 

TTAGGCATCT CCTATGGCAG GAAGAAGCGG 
(2) INFORMATION FOR SEQ ID NO: 2: 

40 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

45 
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(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

5 (vi) ORIGINAL SOURCE: 

(C) INDIVIDUAL ISOLATE: PRIMER ED12, HIV-1-HXB2 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 

10 

AGTGCTTCCT GCTGCTCCCA AGAACCCAAG 
(2) INFORMATION FOR SEQ ID NO: 3: 

15 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 38 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

20 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

25 (vi) ORIGINAL SOURCE: 

(C) INDIVIDUAL ISOLATE: PRIMER ES7 , M13/HIV-1-HXB2 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 

30 

TGTAAAACGA CGGCCAGTCT GTTAAATGGC AGTCTAGC 
(2) INFORMATION FOR SEQ ID NO: 4: 

35 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 39 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



40 



( ii ) MOLECULE TYPE : DNA ( genomic ) 
(iii) HYPOTHETICAL: NO 



45 



(vi) ORIGINAL SOURCE: 
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(C) INDIVIDUAL ISOLATE: PRIMER ES8 r M13/HIV-1-HXB2 
(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4: 

5 

CAGGAAACAG CTATGACCCA CTTCTCCAAT TGTCCCTCA 
(2) INFORMATION FOR SEQ ID NO: 5: 

10 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

15 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

20 (vi) ORIGINAL SOURCE: 

(C) INDIVIDUAL ISOLATE: PRIMER ES 5 , HIV- 1-HXB2 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 

25 

CAATGTATGC CCCTCCCATC 

(2) INFORMATION FOR SEQ ID NO: 6: 

30 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

35 

(ii) MOLECULE TYPE: DNA (genomic) ^ 

(iii) HYPOTHETICAL: NO 

40 (vi) ORIGINAL SOURCE: 

(C) INDIVIDUAL ISOLATE: PRIMER ED4, HIV-1-HXB2 



45 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 
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CACCACTCTT CTCTTTGCCT TGGTGGGTGC 
(2) INFORMATION FOR SEQ ID NO: 7: 

5 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

10 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

15 (vi) ORIGINAL SOURCE: 

(C) INDIVIDUAL ISOLATE: PRIMER ED14, HIV-1-HXB2 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 7: 

20 

TCTTGCCTGG AGCTGCTTGA TGCCCCAGAC 
(2) INFORMATION FOR SEQ ID NO: 8: 

25 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 26 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

30 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

35 (vi) ORIGINAL SOURCE : 

(C) INDIVIDUAL ISOLATE: PRIMER ED5, HIV-1-HXB2 



40 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 8: 

1 

ATGGGATCAA AGCCTAAAGC CATGTG 
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IT IS CLAIMED: 



1. A method of evaluating sequence diversity 
in a mixture of nucleic acids containing a target 
5 sequence, comprising 

selecting amplification primers complementary 
to nucleic acid sequences flanking the target 
sequence, 

combining the nucleic acids and the primers 
10 under conditions that promote the hybridization of 
the primers to the nucleic acids, thus generating 
primer /nucleic acid complexes, 

converting the primer/nucleic acid complexes 
to double-strand fragments in the presence of a 
15 suitable polymerase and all four 
deoxyr ibonucleotides , 

amplifying the number of primer-containing 
fragments by successively repeating the steps of 
(i) denaturing the double-strand fragments to 
20 produce single-strand fragments, (ii) hybridizing 
the single strands with the primers to form 
strand/primer complexes, (iii) generating double- 
strand fragments from the strand/primer complexes 
in the presence of DNA polymerase and all four 
25 deoxyribonucleotides, and (iv) repeating steps (i) 
to (iii) until a desired degree of amplification 
has been achieved, 

denaturing and renaturing the amplified 
fragments to form a population of amplified 
30 fragment DNA duplexes, 

separating the duplexes on polyacrylamide 
gels, and 

analyzing the relative migration of the 
duplexes to establish the relative degree of 
35 sequence relatedness in the population of 
amplified fragments. 
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2. The method of claim 1, where said nucleic 
acid is RNA and where in said converting the 
suitable polymerase is reverse transcriptase. 

3. The method of claim 1, where said 
denaturing is thermal denaturing, and said 
generating is carried out using a thermostable DNA 
polymerase . 

4. The method of claim 1, where said 
separating is carried out by poly aery lamide gel 
electrophoresis, and said analyzing involves 
visualization of the amplification products by 
ethidium bromide staining. 



5. The method of claim 1, where said 
separating is carried out by polyacrylamide gel 
electrophoresis, and said analyzing includes 
transfer of the nucleic acids from the gel to a 
20 support membrane, and hybridization of the nucleic 
acid transferred to the membrane with a labelled 
probe specific for the desired amplification 
products. 

25 6. The method of claim 1, where the primers 

contain at least one detection moiety. 

7/ The method of claim 6, where said 
detection moiety is a radioactive moiety or 
30 biotin. 

8. The method of claim 1, where said nucleic 
acid containing a target sequence is a nucleic 
acid from a microorganism. 



9. The method of claim 8, where said 
microorganism is Human Immunodeficiency Virus 1. 
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10. The method of claim 8, where said method 
Is used to evaluate sequence diversity over time 
in mixtures of nucleic acids containing a target 
sequence, where the mixtures are serially obtained 

5 from a single source. 

11. The method of claim 10, where said 
microorganism is Human Immunodeficiency Virus 1 
(HIV-1) , and where said single source is a patient 

10 infected with HIV-1. 

12. The method of claim 1, where said 
nucleic acid containing a target sequence is a 
nucleic acid from an oncogene. 

15 

13. A method of any of claims 1 to 12, for 
evaluating, in mixtures of nucleic acids, the 
effect over time of a disease treatment, on DNA 
sequence variation of a nucleic acid target 

20 sequence associated with the disease, 

where mixtures of nucleic acids are serially 
obtained from a single source, 

and where the method further includes, 
evaluating the effect of the treatment by 
25 comparing the relative degree of sequence 

relatedness of amplified fragments in each serial 
sample. 



30 



14. A method of evaluating sequence 
diversity between two different sample mixtures of 
nucleic acids, where said nucleic acids contain a 
target sequence, comprising 

selecting amplification primers complementary 
to nucleic acid sequences flanking the target 
35 sequence, 

combining each nucleic acid sample 
individually with the primers under conditions 
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that promote the hybridization of the primers to 
the nucleic acid, thus generating primer/nucleic 
acid complexes, 

converting the primer/nucleic acid complexes 
5 in each sample to double-strand fragments in the 
presence of a suitable polymerase and all four 
deoxyribonucleotides, 

amplifying the number of primer-containing 
fragments by successively repeating the steps of 
10 (i) denaturing the double-strand fragments to 

produce single-strand fragments, (ii) hybridizing 
the single strands with the primers to form 
strand/primer complexes, (iii) generating double- 
strand fragments from the strand/primer complexes 
15 in the presence of DNA polymerase and all four 

deoxyribonucleotides, and (iv) repeating steps (i) 
to (iii) until a desired degree of amplification 
has been achieved, 

mixing together the amplified fragments from 
20 each sample, 

denaturing and renaturing the amplified 
fragments to form a population of amplified 
fragment DNA duplexes, 

separating the duplexes on polyacrylamide 
25 gels, and 

analyzing the relative migration of the 
duplexes to establish the relative degree of 
sequence relatedness among the amplified fragments 
of the population. 

30 

15. The method of claim 14, where said 
nucleic acid containing a target sequence is a 
nucleic acid from a microorganism. 

35 16. The method of claim 15, where said the 

samples of said nucleic acid are collected from a 
number of different geographic locations. 
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17. The method of claim 15, where said 
microorganism is Human Immunodeficiency! Virus l. 

18. The method of claim 14, where the 

5 amplified fragments of one sample nucleic acid are 
labelled with a detection moiety and where said 
labeled fragments are mixed with a molar excess of 
the amplified fragments of the other sample. 

10 19 • A method for detecting the presence of a 

selected nucleic acid target sequence in a nucleic 
acid sample, comprising 

selecting (i) a duplex DNA probe having two 
complementary strands, where the duplex is 

15 homologous to the target sequence and each strand 
contains a detection moiety, and (ii) 
amplification primers complementary to nucleic 
acid sequences flanking the target sequence of the 
nucleic acid, 

20 combining the nucleic acid sample with the 

primers under conditions that promote the 
hybridization of the primers to the nucleic acid, 
thus generating primer/nucleic acid complexes, 

converting the primer/nucleic acid complexes 

25 in to double-^strand fragments in the presence of a 
suitable polymerase and all four 
deoxyr ibonucleot ides , 

amplifying the number of primer-rcontaining 
fragments by successively repeating the steps of 

30 (i) denaturing the double-strand fragments to 

produce single-strand fragments, (ii) hybridizing 
the single strands with the primers to form 
strand/primer complexes, (iii) generating double- 
strand fragments from the strand/primer complexes 

35 in the presence of DNA polymerase and all four 

deoxyr ibonucleot ides, and (iv) repeating steps (i) 
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to (iii) until a desired degree of amplification 
has been achieved, 

mixing the duplex probe with a molar excess 
of the amplified fragments, and 
5 denaturing and renaturing the mix of probe 

and amplified fragments to form a population of 
DNA duplexes, 

separating the population of DNA duplexes and 
the duplex DNA probe on a polyacrylamide gel, and 
10 analyzing the migrations of duplexes^ which 

contain a strand of the probe, relative to the 
migration of the probe duplex, to establish the 
relative degree of sequence relatedness between 
the probe and sample target sequences, 

15 

20 . The method of claim 19, where said 
mixing is carried out at a ratio of 100:1, 
amplified fragments to probe. 

20 21. The method of claim 19, where said 

detection moiety is biotin or a radioactive 
moiety. 

22- The method of claim 19, where said 
25 nucleic acid containing a target sequence is a 
nucleic acid from a microorganism. 

23. The method of claim 22/ where said 
microorganism is a Human Immunodeficiency Virus 1. 

30 

24. The method of claim 23, where said probe 
is selected from the genome of Human 
Immunodeficiency Virus 1. 

35 25. The method of claim 19, where said 

nucleic acid containing a target sequence is a 
nucleic acid from an oncogene. 
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26. A method of any of claims 19 to 25, for 
evaluating the effect of a disease treatment 
procedure on the presence of a selected nucleic 
acid target sequence in a nucleic acid sample, 

where samples are obtained from a single 
source and a first sample is obtained before 
treatment and a second sample is obtained after 
treatment, 

}^-/, a ^^^?^;- ^he-method • further includes, 

^ ^yai^tfn effect of a disease treatment 

procedure on ; tfie pr of a selected nucleic 

ac id^ targe t -s^ufence by coimpar i ng the relative 
degree of sequence^^ the probe 

arid vtarget seggences - of ttie seconci nucleic acid 
sample relative to the Hr nucleic acid sample. 
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